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ABSTRACT 


Brief  consideration  is  given  to  the  possibility  of  using 
alternating  current  asynchronous  interconnections  of  alternating 
current  systems  and  the  desired  properties  of  such  asynchronous 
connectors  are  listed.  It  is  seen  that  a  polyphase  wound  rotor 
induction  machine,  operating  in  the  doubly  excited  mode,  should 
provide  the  desired  properties. 

Differential  equations  are  derived  to  predict  the 
behaviour  of  a  standard  wound  rotor  induction  machine  when  oper¬ 
ating  in  the  doubly  excited  mode.  From  the  differential  equations 
a  program  for  an  electronic  stored  program  digital  computer  was 
written. 

The  transient  currents,  rotor  accelerations,  rotor 
velocities  and  rotor  displacements  are  considered  from  the  moment 
two  systems,  at  slightly  different  frequencies,  are  paralleled 
through  the  machine.  The  same  phenomena  are  studied  when  both 
systems  experience  a  line-to-line  fault,  at  the  machine,  at  the 
moment  of  paralleling.  In  both  cases  it  was  assumed  that  stator 
currents  had  achieved  steady  state  before  paralleling  the  rotor 
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M  Coefficient  of  Mutual  Inductance  (Henries) 

j  A  Complex  Operator  Defined  as  ST 

If  A  Complex  Operator  Defined  as  cos  (***)  4  j  sin(i^) 

0  The  Time  Derivative  =•  ~ 

at 

T  Resistance  (Ohms) 

(f>  The  Electric  Angular  Displacement  Measured  from  the 
Stator  to  Rotor,  Phase  "a"  of  Stator  to  Phase  "a"  of 
Rotor  (Radians) 

W  Work  or  Energy  (Joules) 
p  Number  of  Poles 
J  Polar  Moment  of  Inertia 
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INTRODUCTION 


Two  schemes  are  in  use  for  interconnections  of  synchronous 
alternating  current  systems.  These  schemes  use  either  extra  high 
voltage  direct  current  or  extra  high  voltage  synchronous  alternating 
current.  This  paper  is  concerned  with  a  device  which,  theoretically, 
could  be  a  principal  component  of  an  asynchronous  alternating  current 
interconnection . 

Two  synchronous  alternating  current  systems  would  be 
asynchronously  interconnected  if  each  system  could  be  operated  at  a 
frequency  independent  of  the  other  system,  and  if  power  could  be 
transferred  from  one  system  to  the  other  after  such  an  interconnection 
has  been  made. 

One  possible  type  of  asynchronous  interconnection  is 
composed  of  two  synchronous  machines  coupled  through  a  fluid  drive. 

Such  a  system  is  truly  an  asynchronous  connector  because  the  fluid 
drive  permits  one  shaft  to  rotate  independently  of  the  other  shaft; 
because  the  fluid  drive  permits  transfer  of  mechanical  power;  and 
finally  because  either  synchronous  machine  can  behave  either  as  a 
synchronous  motor  or  as  a  synchronous  alternator. 

The  asynchronous  interconnection  composed  of  two  synchronous 
machines  coupled  through  a  fluid  drive  has  the  following  disadvantages: 

(a)  There  are  two  machines  each  having  nearly  the  same 


iron  and  copper  losses 
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(b)  There  are  two  direct  current  fields  which  must  be  excited 
in  order  to  maintain  synchronism. 

(c)  There  are  two  shafts  rotating  at  synchronous  speed 
and  the  friction  and  windage  losses  are  significant. 

(d)  Fluid  drives  generally  do  not  allow  the  "load"  sides  of 
the  drives  to  rotate  faster  than  the  "motor”  sides  of  the  drives  and 
they  have  not  been  designed  to  transmit  more  than  a  few  megawatts. 

(e)  Synchronous  machines  sustain  fault  KVA  and  any  such 
fault  would  be  supplied  by  two  electric  systems. 

The  writer  sought  a  device  which  would  have  the  properties 
of  an  asynchronous  connector  without  the  above  disadvantages .  The 
ideal  alternating  current  asynchronous  connector  would  have  the 
following  properties : 

(a)  One  shaft  would  revolve  at  slip  speed. 

(b )  The  device  would  not  convert  electric  power  to  mechanical 
power  at  any  point  in  the  connector. 

( c )  The  asynchronous  device  would  not  require  direct  current 

excitation. 

(d)  A  fault  on  either  electric  system  would  be  electrically 
isolated  from  the  other  system  by  the  asynchronous  connector. 

In  the  conventional  wound  rotor  induction  motor  there  are 
four  principle  components: 

(a)  A  stationary  part  called  the  stator  has  nearly  identical 
phase  resistances  and  has  nearly  identical  phase  reactances.  The  stator 
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phasa  windings  are  nearly  symmetrically  placed. 

(b)  A  rotating  part  called  the  rotor  has  nearly  identical 
phase  resistances  and  has  nearly  identical  phase  reactances.  The 
rotor  phase  windings  are  nearly  symmetrically  placed. 

(c)  A  set  of  slip  rings,  to  which  the  rotor  windings  are 
connected,  and  a  set  of  brushes,  to  which  the  slip  rings  electrically 
connect,  are  used  to  insert  impedances  into  the  rotor  circuits. 

(d)  External  impedances  of  nearly  equal  phase  values  are 
inserted  into  the  rotor  circuit  to  control  rotor  speed. 

Doubly  excited  wound  rotor  induction  machine  operation  is 
said  to  exist  when  excitation  from  separate  sources  is  applied  to  the 
rotor  and  stator. 

If  the  rotor  of  a  doubly  excited  wound  rotor  induction 
machine  rotates  at  the  slip  speed  the  rotor  fluxes  are  in  time  phase 
with  the  stator  fluxes  and  the  mutual  fluxes.  Under  such  circumstances 
the  rotor  speed  will  be  given  by  the  formula 
ns  - 

where  is  the  rotor  speed  in  revolutions  per  minute 

p  is  the  number  of  poles  in  both  stator  and  rotor 

f,  is  the  frequency  in  cycles  per  second  applied  to  the  stator 

$2  is  the  frequency  in  cycles  per  second  applied  to  the  rotor. 

Clearly  -f,  and  can  be  different  and  independent.  Therefore 
the  first  criterion  of  an  asynchronous  connector  is  satisfied.  Further, 


t 

<  t  .  i„  ■' 

•  c  '  :  j  :  ;o 


...  •;  .  j  . 

■ 


:  r3  , 

* 

«  ■  f  '  ■  1  ■  ,  | 


IV 


if  the  rotor  rotates  at  slip  speed  and  the  ratio  of  stator  to  rotor 
voltages  changes,  the  air  gap  flux  density  changes  and  a  change  in 
magnitude  of  KVA  flow  results.  Generally  speaking  this  implies  a 
change  in  power  flow  and  a  change  in  KVAR  flow. 

The  objects  of  this  research  are  as  follows; 

(a)  The  steady  state  advantages  and  disadvantages  of  a 
doubly-excited  wound  rotor  induction  machine  are  fairly  easily 
studied,  but  the  starting  phenomena  are  net  so  readily  studied 
and  the  starting  phenomena  form  the  main  subject  of  this  research. 

(b)  The  fundamental  problem  is  to  select  a  suitable 
asynchronous  connector.  This  research  indicates  that  the  conventional 
doubly-excited  induction  machine  is  likely  to  be  unstable. 

(c)  The  advantages  and  disadvantages  of  asynchronous 
interconnections  are  numerous  but  difficulties  can  be  resolved  only 
through  increased  interest  in  the  subject.  One  object  of  this  paper 

is  to  stimulate  enough  interest  that  the  disadvantages  may  be  overcome. 
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CHATTER  I 


KIRCHHOFF'S  LAW  AND  DERIVATION  OF  VOLTAGE  EQUATIONS 

I.  INTRODUCTION 

The  following  considerations  are  necessary  before  analysis 
can  be  attempted: 

(a)  One  mast  be  able  to  determine  the  instantaneous  currents 
and  voltages  in  both  the  stator  and  rotor. 

(b)  One  must  be  able  to  predict  the  position  of  the  rotor 
and  its  velocity  at  any  instant  in  order  to  determine  the  electrical 
torque  induced  and  the  retarding  mechanical  torques. 

In  this  chapter  derivations  are  presented  which  enable  one 
to  predict  the  currents  flowing  at  any  instant.  Induced  torques  are 
studied  in  Chapter  II. 


Throughout  this  paper  the  machine  is  assumed  to  be  connected 
as  in  figure  1.  The  rotor  body  is  assumed  to  rotate  at  angular  speed 
electrical  radians  per  second  and  positive  a)  is  assumed  to  be 
counter-clockwise . 
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II.  RESTKICTICNS 

All  of  the  formulae  in  Chapters  I  and  II  are  dependent  on 
the  assumptions  that: 

(a)  Stator  and  rotor  have  negligible  hysteresis  and  eddy 
current  losses. 

(b)  The  air  gap  is  uniform  and  the  air  gap  flux  is 
sinusoidally  distributed  in  space. 

(c)  The  stator  and  rotor  are  perfectly  symmetrical  with 
respect  to  resistance  and  inductance  as  well  as  physical  location 
of  the  phase  windings. 

(d)  The  resistances  and  inductances  of  stator  and  rotor 
windings  are  independent  of  current  magnitude  and  frequency. 

(e)  The  se If -inductances  of  stator  and  rotor  windings  are 
independent  of  their  relative  positions  whereas  the  mutual  inductance 
between  stator  winding  and  a  rotor  winding  varies  as  the  cosine  of 
the  electrical  angle  between  them. 

(f)  The  capacitances  produced  by  proximal  locations  of 
windings  one  to  another  have  negligible  effects. 

(g)  The  impedance  of  sources  of  potential  for  stator  and 
rotor  are  negligible  and  hence  the  stator  and  rotor  may  be  regarded 
as  being  connected  to  semi-infinite  sources. 
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III.  THREE  PHASE  WINDINGS 


The  "a"  phase  flux  linkage  which  is  the  result  of  currents 
in  phases  "a",  "b"  and  "c"  under  the  assumptions  in  Section  II  is 
given  by 

=  ‘-a  V  *  Mab  ‘b  *  Mac  <’> 

where  L  ^  is  the  self-inductance  of  phase  "a" 

/V\ab  is  the  mutual  inductance  between  phase  "a"  and  phase  'b" 
Mac  is  the  mutual  inductance  between  phase  "a"  and  phase  "c" 

and  where 


inductances  are  in  henries 


1 


currents  are  in  international  amperes  . 


With  the  assumptions  of  Section  II  being  made, 

A 

r  Mbc. 

By  Kirchhoff's  current  law,  flux  through  phase  "a"  caused  by 
positive  current  in  phase  "b"  would  be  oppositely  directed  to  flux 
through  phase  "a"  caused  by  positive  current  in  phase  "a"  itself. 

Also  by  Kirchhoff's  current  law 

‘‘b  +  6C  =  ~  CCL 

and  in  consequence  of  the  formula 


1  Throughout  this  paper  the  MKS  system  of  units  is  used. 
If  CGS  units  were  used,  the  units  would  be  abhenry  and  abampere 
respectively.  The  potential  would  be  in  abvolts  for  CGS  units  and 
volts  for  the  MKS  system. 
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it  is  obvious  that 

=  (  i<L  -  Mafc>  )  4 «.  (2) 

Equation  (2)  is  valid  provided  no  "neutral1'  current  flows. 

If  neutral  currents  flow,  equation  (1)  becomes 

^ a,o  ~  L'a*  ^a*o  ^  ^  o-b  ^  be>  ^ac  ^  c  o  (3 ) 

Since  ibo  -  iCe>  *  Ca.0  by  definition 

^  CL  0  ~  (  CL.  +2  M  a.  b  )  t*  4,  a  (a  cl) 

Since  the  flux  through  phase  "a"  caused  by  positive  current  in  phase 
"b"  is  oppositely  directed  to  flux  through  phase  "a"  caused  by  positive 
current  in  phase  "a"  itself,  and  Mao  in  equation  (1)  must  be 

assumed  to  be  negative.  If  and  Ma<.  ere  assumed  negative  then 
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IV.  MUTUAL  FLUX  LINKAGES 


For  three  phase  windings  the  mutual  flux  linkage  with  rotor 
phase  "a"  caused  by  currents  in  a  symmetrical  three  phase  stator 
winding  are  as  follows; 

^(3cl  =  ^ p&da,  ^</ol  4  MjSaocb  Lo<b  4*  ^poLdc  C4) 

Equation  (4)  is  valid  for  all  degrees  of  balanced  or  unbalanced  three 


phase  circuits  . 
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Assuming  sinusoidal  flux  distribution  in  the  air  gap,  the 

following  equations  are  obvious: 

, .  <  .  i  j  _  j  d , 

Mfra.da,  -  M  cos  <f>  -  M  (  Q  +  e  ) 

2 

M  t>  -  M  COS  (<£  -  2J1)  =  ^  (  e  ^  e  J 

3  2 

»  ./  /A  -  i  /A 

(see  ■£;%.  Z  ) 


-  m'  (^2e^y.  ye'J^) 
z 


simi larly 


MpaoCc  =  —  (  yr  t  y2  e~J*> ) 

2  7T/3  2  7773  ^ 

- - >■ 


(  See  /,<».  2 ) 
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by  definition 


lo.a  =  lcu  +  L<xz  +  L  (to 

bdb  -  Y  Ldt  +  ¥Ldt  *  Lo io 

4 tic  -  VCQCI  *  Y  C&Z  +  LUO 

Substitution  into  equation  (4)  yields 

^jja  =  ~  (  e  +  e  ^  ^oLi  +  t'dZ  + 

+  £  (  T  e  t  ye  )\  y  lu,  +  yctiz  +  i«0) 

2 

+  (A  (ye^  +  y%e~J*)(  Yc<t,  +  Yxi^z  t  i*c) 

z 

This  may  be  shown  to  be  equivalent  to  the  following  equation: 

Vo.  s  J-  m'  (  i*.  e~J *  +  lv2  } 


{S) 


Similarly  the  following  equation  may  be  deduced: 

_  -  j<p 

m  (t4l  e 

2 


^  2.  m'  ( e  '  v-  LpZ  e 


) 


(5a) 


-  .  i.  )  )  C  i 


6. 


Substituting  M  =  1m  into  equation  (5)  and  (5A) 

2 


A\  (  CKI  e  4 

L  V2 

J  ) 

e  > 

(6  ) 

'Wa.  » 

Art,  (  la/  3  -f 

-j<6  s 

e  ) 

(7) 

Thus  if  the  flux  is  sinusoidally  distributed  in  space  the 


mutual  flux  linkage  is  independent  of  the  zero  sequence  currents  in 
either  stator  or  rotor  since  no  zero  sequence  components  occur  in 
equations  (6)  and  (7).  While  no  zero  sequence  components  occur  in 
the  air  gap  it  is  not  obvious  that  zero  sequence  terms  do  not  occur 
in  the  equation  of  torque.  The  study  of  torque  will  show  that  such 
is  true.  Torques  are  studied  in  Chapter  II. 


V.  VOLTAGE  EQUATIONS 


Displacement  currents  and  leakage  currents  through 
insulation  are  small  and  completely  negligible.  Such  being  the 
case  the  fall  in  potential  between  the  terminals  of  a  stator  or 
rotor  consists  of  two  components  as  follows: 

(a)  The  resistance  drop  r 0  dissipates  power  and 
provides  electromagnetic  damping. 

(b)  The  variation  of  total  flux  linkage,  f!_*  ,  generates 

dt 

a  voltage  but  dissipates  no  power. 


and  (7) 


The  total  flux  linkage  consists  of  two  parts: 

(a)  The  self -flux  linkage  is  given  by  equation  (2)  and  (4) 

(b)  The  mutual  flux  linkage  is  given  by  equations  (6) 


'?)  roMfjrp* 

. 

. 

, -<9va  O'  :  otf  '  '  ■  ’-fr  ref  "o' 

tot  °fp'>  lioy"';  r '  -  ■  >f*  ;.I  xuXTt-1  ".os  e iT 


.( T)  baa 
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The  fall  in  potential  between  the  terminals  of  phase  "a" 
of  the  stator  caused  by  positive,  negative  and  zero  sequence  currents 
is  as  follows: 

+  ^Z  +  ^  4  ^  +  +  ^OiO  ) 

4-  AA  D  (  eJ<Z*4  e'j^)  <3) 

The  zero  sequence  flux  linkages  do  not  appear  in  the  above  equation. 

The  corresponding  equation  for  phase  "b"  of  the  stator  is  as  follows: 

^  y  =■  If  *  V  Kz  =  ^  i  <v/  •*  V  taz  +  ) 

4  M  p  (  yj*  f  ipi.  y  e.  J<*j  (9) 

Similarly  for  phase  "c" 

r  V*  ,L  W  °&2.  +  ^<3  =  ■*  tvD)(  ip  Lyt  4  P  yz.  +  ^tyc  ) 

■i  md  C  lpt  +l^  y-'e  J<^)  (lo) 

where  0  radians  has  been  defined  (see  fig-  2) 

0.  *x  radians  is  the  electrical  angle  from  "b*1  phase  of  stator 
to  "a"  phase  of  rotor 

A-  4ir  radians  is  the  electrical  angle  from  "c:i  phase  of  stator 

3 

to  "a"  phase  of  rotor. 

One  observes  that  addition  of  equations  (8),  (9)  and  (10) 
causes  cancellation  of  all  terms  except  those  of  zero- sequence .  When 
this  sum  is  simplified 

The  negative- sequence  components  of  potential  and  current 
can  be  eliminated  by  multiplication  of  equation  (8)  by  p  and  then 
by  subtraction  of  equation  (9)  from  the  result.  Thus: 


.  ..  '•  '• 


8. 


Ud)  tUl  +  02) 

Similarly  the  negative  sequence  equation  is  as  follows: 

V«x  -  (1*  +  L*  D)  La2  +  MD  (4tt  )  (/3) 

The  corresponding  rotor  equations  can  be  deduced  in  a 


similar  manner  and 

they  are  as 

follows : 

+  MD  (  lui  e'J *) 

0  4) 

V  = 

(  +  Lp  D) 

+  MD  c  i„a  e’J*) 

(is) 

and  - 

(  +  Lfie> 

(|6) 

In  general  the  zero  sequence  components  are  negligible 
since  the  earth  leakage  impedance,  of  very  large  magnitude,  is  in 
series  with  the  transient  impedances  above  derived. 

While  of  no  practii  purpose  unless  the  star  point  of 
either  stator  or  rotor  is  earthed,  equation  (16)  was  derived  to 
demonstrate  that  zero  sequence  transient  impedance,  (  +  L^0D  )  , 

is  very  small  and  to  demonstrate  the  necessity  for  not  earthing  the 
star  point  of  either  stator  or  rotor  of  a  doubly  excited  induction 


machine . 


<  t 

< 

t 

<  < 


CHAPTER  II 


TORQUE 

I.  INTRODUCTION 

Whenever  a  current  carrying  conductor  links  a  magnetic  field, 

a  force  is  produced  which  attempts  to  minimize  the  linkage  of  flux. 

Several  elegant  methods  have  been  evolved  which  predict  the  force 
1  2 

produced.  9  Some  methods  employ  tensor  calculus  to  derive  the  force 
produced  and  others  consider  the  torque  by  applying  clever  artifices. 

In  the  end,  all  methods  employ  the  same  assumptions  and  all  employ 
methods  derived  from  energy  considerations.  Since  the  methods  employ 
the  same  assumptions  and  devolve  from  energy  considerations,  little 
can  be  gained  by  employing  tensor  notation  when  more  easily  understood 
methods  can  be  employed.  In  the  bibliography  of  the  books  below  are 
selected  references  for  explanation  of  the  torque  phenomenon  based  on 
tensor  notation.  Unfortunately  none  of  the  references  explain  the 
torque  phenomena  in  terms  of  symmetrical  components.  To  explain  the 
torque  phenomena  in  terms  of  symmetrical  components  requires  considerable 
knowledge  of  symmetrical  spinor -matrix  calculus.  This,  then,  is  the  main 
reason  for  abandoning  tensor  notation. 


1 

L.  V.  Bewley  Tensor  Analysis  of  Electric  Circuits 

and  Machines,  The  Ronald  Press  New  York  (1961) 

2 

G.  Kron  Tensors  for  Circuits,  Dover  Publications,  Inc. 

New  York  19^-2,  copyright  1959  Dover  Publications. 
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II.  TORQUE  FOR  SINGLE  PHASE  MACHINE 


Considering  a  single-phase  rotor  single-phase  stater,  the 
total  stored  electromagnetic  energy  of  such  a  system  is 

w=  ±L*C+  +/v\  joules 

If  one  assumes  that  the  stator  and  rotor  coils  are  connected  to  sources 
of  electromotive  force  which  maintain  the  currents  and  ^  constant 
and  if  the  rotor  moves  a  distance  dx-in  time  dt  ,  the  increase  in 
stored  electromagnetic  energy  is  as  follows 

dw  =  l^dAA  joules 

where  and  Lp  have  the  same  sign  if  they  produce  magnetic  fields  in 
the  same  direction  when  the  displacement  is  zero. 


Electromotive  forces  are  generated  in  all  coils  by  the 


variation  of  the  mutual  flux  linkages  and  are  as  follows: 

e 


&IX  -  —  _d  (  tp/vy  ) 
dt 


=  -  Lr. 


dM 

dt 


-  -  A  (  Lgc  /^) 

dt 

r:  -  C  oc  d  <V| 

dt 


The  combined  electrical  output  of  rotor  and  stator  is  therefore 

(e*  watts  and  in  the  time  dt  the  energy  output  would 

be  ( e*  t’oc  +  tp  )  dt  =  -  2  i u  l ^  d aa  joules.  The  work  done  by  the 

•  , 

applied  force  F  is  F  dx.  joules.  Thus  the  total  work  dene  by  the 

i 

applied  force  F  is  given  by 

f'c/x  =  (increase  in  stored  energy)  +  (energy  delivered) 


< 


- 


• 

l  . ) 

11. 


Newton's  Third  Law  dictates  that  the  applied  force  must  be 
opposed  by  an  equal  developed  force  or  that  the  developed  force,  F, 
is  given  by  the  equation 

F  -  -  F  ’ 

Thus  Pd*  -  joules 

Hence  F  „  dM  newtons 

1  dx 

Since  torque  can  be  determined  by  multiplying  force  by  effective 
radius  then 

T  =  dM  .r)  :  dM  newton-meters. 

I  rd0  J  d<f> 

If  stator  and  rotor  each  has  p  poles  then  the  rate  of  variation  of 

M  with  respect  to  an  angular  displacement  0  must  be  P.  times  as 

2 

fast  as  if  only  2  poles  are  in  the  machine.  Under  such  circumstances 

T  _  I  p  l*  bA  dM  newton-meters 

“2  d  0 

where : 

all  angular  measurements  are  in  electrical  radians 
p  is  the  number  of  poles  on  stator  and  rotor. 

<p  is  as  defined  in  the  table  of  symbols  and  Fig.  2. 

In  the  "ideal"  machine  the  mutual  inductance  between  a  stator 

and  a  rotor  phase  is  given  by  the  equation 

i  .  •  ,  J<t>  -j<£> 

M  _  M  cos  d>  -  J.  M  (  e  +  <2 

z 

Differentiating  with  respect  to  0  yields 

T  -  1  p  Lu  L  j,  M.  (  -  e  ) 

a  1  6  Z 


07) 


.  :  • 

«  <  '  i  j  'j-  .  .. 


•  * 

. 

-  J 

: 
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III.  THREE  PHASE  WINDINGS 


The  three  phase  stator  currents  induce  torques  on  the  phase 
"a"  windings  of  the  rotor.  Since  a  three  phase  stator  would  have 
three  separate  currents,  three  equations  similar  to  equation  (17)  must 
be  combined.  The  torque  on  phase  "a"  windings  of  the  rotor  caused  by 
currents  in  phase  "a"  of  the  stator  is  given  by 

-r  '  J  &  -  J  &>  \ 

Tac  s  }  Z  (  e  -  e  } 


^  I  )  (  (,#,  V-  Lg 2  +  6(X0  )  (  e  -  &  ) 


2  Z 


(n  a) 


The  torque  on  phase  "a"  windings  of  the  rotor  caused  by  currents  in 


phase  "b"  of  the  stator  is  given  by 


T ^  Z  ^  (■  Lft  *  Lfiz  +  ^/So)(  V  t'OU  +  'YC&Z  y  C&o)  ( e 


j'ct-z?)  -jv-m 


-  e 


—  ^  Z.  Ai  (  ^  +  ^52  +  ( Iff  I  Y  L tfz  +  i do  }(  P  C  -  T]/"  G  )  ( I  7  B) 


x  2 


The  torque  on  phase  "a"  windings  of  the  rotor  caused  by  currents  in 
phase  "c"  of  the  stator  is  given  by 

_  .  •  N/  •  i.  .  w  J(t-H) 

'ac  =  y  (  tp,  4  z  Y  £&/  4  ff  £*0(2  +  t'cto*  ^  ^  ^  ' 


-  ^  .E  m  ( tp,  v  +  c^oY  Y toe/  +  Y  *  ^«o)  (  y Q  -  p/  e  ^ )  07 c) 

The  torque  on  phase  "b"  windings  of  the  rotor  caused  by  currents  in 


phase  "a"  of  the  stator  is  given  by 


Tba=  !*  (.ir'Ljn  +  Vifiz+ir 


f  +  7/f  ^  ^  OL  /  ^  ^(X  2  71  ^cfo)  ^  ' 


=  j  E.  ^  ^  <^2  ■’t  t'/BoY  +  L-oiz  +  ^ao)  (  a  -  e  )  ( IJD) 


■  j<p 


z  2 


The  torque  on  phase  "b"  windings  of  the  rotor  caused  by  currents  in 


. 

. 
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phase  "b"  of  the  stator  is  given  by 

jV 


Tt»b  =.  ^  Z  HQ.  ‘-c/b  ^ 


e  -  e 


j  0 


=  i  £  %  (Vz^,  +  &/*>)(  v\,*r*<rz  +£«c)(eJ<*-  (I7E) 

The  torque  on  phase  "b"  windings  of  the  rotor  caused  by  currents  in 


phase  "c"  of  the  stator  is  given  by 

Tbc=  £  (V  +  +  +  0(e  ~e  ) 


=  \  £.  £  ( y'LA,+  yU*  y  +  y\t  +  ^vc)  (  v'eJ(*>  -  pre  *'*') 

'2  2.  ~  '  ' 


'}<*> 


(>7F) 


Similarly  one  can  deduce  three  more  equations  for  the  torques  on 
windings  of  phase  "c"  of  the  rotor  caused  by  the  currents  in  phases 
a,b,c  of  the  stator.  These  equations  are  as  follows: 


Tca=  i  £  £  ipc  i*a.  <  e 


) 


2  Z 


.  «  «  2,  .  p  J  (£>  J  0  \ 

z  h  —  —  (V Lp,  iff  632  +  i/s«)(  Lon  •*  L<*2.  +  ^#0)  (  V  ®  -  y/  g  ) 

,  •  .  j +  -i(4  +  n 

Tcb-.tP^  *  -  e  ^  3  > 

*  2  2  P 


=  a  p  aa'  ( ?/**,  *  *  t/io)(yr2i<x,  +  ty igi  +  4oX  f/e  -  v'e~J*) 

*  2  2 

T.  _  «  /  -J0\ 

n ^  —  £0.  pc  ^«c  (  ®  ®  ' 


cc 


2  2 


=  j.|  ^  +V2L/>2  +  V*i«z  +loto)(.eJ*~  e~J<6) 


(17^) 


0?h) 


(I'M) 


Addition  of  equations  ( 17A) ,  ( 17B) ,  (17C)5  ( 17D) ,  ( 17E) ,  ( 17F) ,  ( 17G) , 
( 17H)  and  (171)  gives 


T  =  |  *  pM  (  e  -  i«,  ifit  e  ) 


(is) 


The  zero  sequence  currents  have  no  direct  effect  on  the  torque 


developed  since  zero  sequence  currents  do  not  occur  in  equation  (18). 


'  '  «  '  *  <  c  •  ••  )  «  <  i  i 


14. 


This  is  reasonable  since  no  flux  exists  in  the  air  gap  and  zero 
sequence  currents  in  the  rotor  are  in  time  phase  and  hence  can 
produce  no  torque. 

With  regard  to  equation  (18)  the  following  points  are 

observed : 

(a)  the  sign  of  the  exponent  for  the  exponential  function 
is  opposed  to  the  sequence  of  the  stator  current  which  multiplies 
the  exponential  function. 

(b)  the  sign  of  the  negative  sequence  stator  current  is 
positive  while  that  of  the  positive  sequence  current  is  negative. 

(c)  zero  sequence  currents  do  not  appear  in  the  torque 

equation. 


The  above  derivation  is  a  very  slight  modification  of  the 
derivation  as  presented  by  Professor  W.  V.  Lyon  in  section  3.6  of 
his  text,  Transient  Analysis  of  Alternating  Current  Machinery. * 

The  modification  is  the  inclusion  oi  zero  sequence  currents.  The 
fact  that  zero  sequence  currents  were  not  explicitly  included  by 
Professor  Lyon  is  probably  because  they  do  not  occur  in  the  torque 
equation.  This  writer  includes  them  in  order  to  clarify  his  concept 
of  the  torque  phenomenon. 

Since  zero  sequence  effects  do  not  have  any  bearing  on  the 


1 


W.  V.  Lyon 


Transient  Analysis  of  Alternating  Current  Machinery 
The  Technology  Press  of  M.I.T. 

Wiley  1954. 
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instantaneous  torque  developed  and  since  the  earth  leakage  impedance 
of  a  wound  rotor  induction  machine  is  very  large  one  can  conclude 
that  the  zero  sequence  currents  would  be  very  small. 

For  these  reasons  zero  sequence  components  are  neglected 
henceforth. 


. 


. 


CHAPTER  III 


THE  MECHANICAL  EQUATIONS  OF  MOTION 
The  mechanical  equation  of  motion  is  given  by  the  differential 


equation : 


J  d(A  4.  R, 60  =  ta 

dt 

2  - 1 

where  J  is  the  polar  moment  of  inertia  (ft.  lbs.  sec  radian  ) 

w 

_2 

is  the  rate  of  change  of  angular  velocity  (radian  sec  ) 

dt 

R  is  the  damping  function  (ft.  lbs.  sec  radian  ) 

w 

4)  is  the  angular  velocity  (radian  sec"1) 

Ta  is  the  net  applied  torque  after  "static"  or  kinetic 

friction  has  been  deducted  (ft.  lbs.  ) 

w 

If  t  is  used  to  designate  the  interval  of  time  under  consideration 
and  i-t  is  used  to  designate  the  time  interval  previously  considered, 
then  equation  (1)  can  be  expressed  as 

±  (  Ta.  -  Rwp() 

dt  J 

Multiplying  equation  ( 1A)  by  a  small  time  increment,  At  , 
gives 

dt 


(0 


OA) 


d6)£  At  ~  A  CO  -  ~  AC  (  TA.  - 

J 


Then 


s'  <Aa  t  E  a -  <tt0  ■*  Z  (  ta.  -  R  &*/.,) 
i -o  i  c o  J 

Multiplying  equation  (2)  by  a  small  time  increment,  At  , 


At 


(IB) 


(2) 


gives 

Cdt  At  c?  A  i;  'l 

n 

(  ti0  +  Z 

Then 

n 

L  1*0 

n 

S,  A  Si 

l  ZD 

=  t  {A 

i- o 

(3) 
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Assuming  that  the  torque  can  be  calculated,  and  assuming 
that  R  and  J  can  be  determined,  equation  (2)  predicts  the  rotor 
velocity  while  equation  (3)  predicts  the  rotor  displacement. 

The  rotor  velocity  is  necessary  in  determining  the  induced 
voltage  drops  and  the  rotor  displacement  is  necessary  in  determining 
the  torque. 

With  reference  to  fig.  1  below,  one  can  easily  deduce  the 


value  of  0  in  electrical  radians. 


ai  t*  o  at  t-  1 1  at  t  -  1 t 

F,q.  /  Flux  Displacements  u>ih  Locked  Rotor 

displacement  between  the  instantaneous 
the  instantaneous  rotor  flux  phase  "a" 

Ojp  6  -  £  ~  ^oc7-  ~~  S  o 

or  (cdp  -  cou )  t  +  80  -  S 

For  obvious  reasons,  then,  one  deduces 

$  ~  S  ~  (  F  S  <5 


If  S0  is  the  physical 
displacement  between  the 
instantaneous  stator  flux 
phase  "a"  and  the  instantaneous 
rotor  flux  phase  "a"  at  t=0, 

Met*  CDpt  -  27 r/pt 

and  £  is  the  physical 
stator  flux  phase  "a"  and 
at  time  t=t,  then  fiom  Fig.  1 

that 


provided  the  rotor  body  does  not  rotate.  In  general,  however,  the 
rotor  body  rotates  and  consideration  must  be  given  to  this  fact.  With 
reference  to  fig.  2  one  uses  the  following  notation: 


e  I  .  •  >  .  •  i  >'•  i  iy.  : 

<  .  ,  1  ,  . 

. 

.  .  .  ■ . 


. 


< 


•  •  .  .  t 


e 


<  c 

«  < 
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S d  is  the  "physical11  displacement  between  instantaneous 
stator  flux  "a"  and  the  instantaneous  rotor  flux  phase  "a"  at  t=0 
6 .)**=■  }  cJpt  »  2*fpt 

is  the  "electrical"  displacement  between  instantaneous 
stator  flux  "a"  and  the  instantaneous  rotor  flux  phase  "a"  at  t=t 
and  S  is  the  displacement  of  the  rotor  phase  "a"  in  time  t. 


6  =  0  6  =  0  t=  t,  at  t  =  t,  at  t=-t. 

Fig.  2  F/ok  Disg!  a.ce  merits  wit/?  Unlocked  Rotor 


From  Fig.  2  one  can  deduce  that 

Sa  +  S  4  Cl)  p  £  ~  $  -  a)#  t 

or 

(cJp-  to*)  t  t  (S0 1  $)  =  6 

Clearly  SQ  t  S  can  be  replaced  by  S  if  Sa  is  assumed  as 
a  constant  of  integration.  Then 

<f>  =  ztt  (  fp  -  4,  )  t  +  S  (4) 

Differentiating  equation  (4)  with  respect  to  time  gives 

<p  _  djfi  _  z*  dj> 

dt  d  t 


(5) 
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Thus : 
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CHAPTER  IV 


DAMPING  FUNCTIONS  AND  POLAR  MOMENTS  OF  INERTIA 

I.  DAMPING  FUNCTIONS 

Ideally  friction  and  windage  would  consume  no  power  and 
require  no  torque.  In  actual  fact,  friction  and  windage  give  rise 
to  damping  functions.  Since  these  damping  functions  are  important 
to  the  whole  study,  the  damping  functions  are  studied. 

Frictional  forces  do  not  occur  until  motion  is  on  the 
verge  of  occurring  or  actually  occurs.  At  such  time,  as  motion 
occurs  the  frictional  forces  suddenly  occur  and  these  frictional 
forces  are  finite.  Clearly,  however,  the  rate  of  change  of  the 
forces  with  respect  to  time  is  not  defined,  and  neither  is  the 
rate  of  change  of  the  forces  with  respect  to  speed. 

Since  no  power  is  consumed  until  motion  actually  occurs 
one  can  assume  the  following: 

(a)  After  motion  actually  occurs,  the  component  of  friction 
windage  and  viscosity  power  loss  that  is  linearly  dependent  on  speed 
must  be  purely  frictional. 

(b)  After  motion  actually  occurs,  the  component  of  friction 
windage  and  viscosity  power  loss  that  is  dependent  on  the  square  of 
speed  should  be  composed  of  windage  and  viscosity. 

Clearly  the  total  power  consumed  in  friction,  windage  and  viscosity 
should  be  continuous  except  at  the  origin  (0  speed,  0  power) ,  and 
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ail  of  its  derivatives  should  be  continuous. 


Mathematicians  have  shown  that  any  continuous  function  can 
be  expanded  in  a  power  series.  Such  proofs  are  found  in  most  textbooks 
of  advanced  calculus.  If  the  assumption  is  made  that  friction  end 
windage  power  is  continuous  with  respect  to  speed  then 


n~  oo 


^F-tW  - 


Pn  5 


0) 


However 
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r F  ■+  uo 


-  c 


cut  s  -  O 


(IA) 


Consequently 

=0 

Hence  00 

=  zr  Pn  s  ” 

n  - 1 

Ordinarily  the  series  in  equation  ( 1A)  converge  very  quickly 

3 

and  one  can  forget  terns  such  as  P  5  . 


Referring  to  Equation  (1)  Chapter  III, Rea  must  have  the 
dimensions  of  torque  and  hence  R  must  have  the  dimensions  of  torque 
seconds  per  radian.  Equation  ( 1A)  above  can  be  modified  as  follows: 

CO 

■*  w.  -  Tf+w  -  Z]  Pn  5  (2) 

5  n  -  * 

Also  Rttcan  be  modified  as 

P^d  _  R.  2  F'S 

<oO 

From  (2)  one  can  readily  deduce  that 

ETTS  R  _  R,  S 

GO 

or 

R  =  60  (3) 
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If  the  friction  and  windage  losses  can  be  determined  at  a 
reasonable  number  of  speeds  and  if  a  quadratic  regression  can  be 
applied,  then  R  can  be  determined  from  such  a  technique.  The  writer 
attempted  such  a  series  of  measurements  and  the  technique  employed 
is  outlined  in  Section  III  of  this  paper. 


II.  POLAR  MOMENTS  OF  INERTIA 


The  angular  kinetic  energy  work  theorem  states  that  the 
kinetic  energy  possessed  by  a  body  rotating  at  co  radians  per  second 
is  i  Jo)2-  H  -  ibs^  ,  where  J  is  the  polar  moment  of  inertia.  If  all 
supply  of  power  to  such  a  rotating  device  is  removed,  then  the  rotor 
will  be  retarded  by  friction  windage  and  viscosity. 

Hence 
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Presumably  equation  ( 1A)  Section  I  has  enabled  one  to  plot  the  locus 
of  friction,  windage  and  viscosity  power.  Then  the  area  between  the 
axis  of  speed  and  the  locus  of  friction,  windage  and  viscosity  power 

rt  n=°° 

is  /  Z7  Pn  s  as  .  This  last  statement  is  intended  to  demonstrate 
Jo  o-i 

the  complexity  of  the  problem  since  the  increments  of  integration  do 

not  agree  and  since  5  =  d  &  f  d 

eft  \dt 


Unfortunately  5  is  a  function  of  time,  and  such  terms  as 
/  cannot  be  defined  analytically.  Neglecting  such  terms 
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would  be  tantamount  to  assuming  that  no  damping  torques  exist.  Clearly 
other  approaches  are  necessary. 


III.  RETARDATION  TECHNIQUE  FOR  DETERMINING  POLAR  MOMENTS  OF  INERTIA 

If  one  applies  the  differential  equation 

Jo)  +  R6)  =  O  (4) 

one,  in  effect,  assumes  that  no  purely  frictional  torques  exist  since 
R  is  windage  and  viscosity  dependent.  Appendix  C  contains  several 
figures  which  indicate  that  purely  frictional  torques  do  exist.  The 
retardation  technique,  which  is  the  standard  technique,  assumes 
differential  equation  (4)  to  be  correct. 

If  the  body  is  rotating  at  rate  radians  per  second  when 
all  supply  of  power  fails  then  the  solution  to  equation  (4)  is  as 
follows 

t/j 

cd  _  40o  e  (s) 

If  more  moment  of  inertia  is  added  to  the  revolving  body 
the  following  complications  arise: 

(a)  The  frictional  force  in  the  bearings  increases  because 
of  the  increased  bearing  pressure. 

(b)  The  windage  force  increases  because  an  increased  area 
is  added  to  come  in  contact  with  air. 

(c)  The  viscosity  effects  in  the  bearings  are  unknown  for 


the  following  reasons: 
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(1)  The  oil  films  become  smaller  and  shearing  forces 
in  the  oil  films  become  more  dependent  on  actual 
bearing  temperature. 

(2)  The  thickness  of  the  oil  films  govern,  to  some 
extent,  the  coefficient  of  viscosity  of  the  oil. 

(3)  With  very  thin  oil  films  the  probability  of  metal 
to  metal  contact  is  increased. 

(d)  The  center  of  gravity  of  the  shaft  moves  axially  and 
results  in  non-equal  bearing  pressures  and  consequent  increase  in  the 
probability  of  metal  to  metal  contact. 

Clearly  some  dubious  assumptions  are  involved  in  use  of  the 
Retardation  Technique  since  the  effects  (a),  (b),  (c)  and  (d)  are 
neglected.  Assuming,  also,  that  the  purely  frictional  effects  are 
negligible  one  can  repeat  the  experiment  with  a  moment  of  inertia 
increased  to  j'  .  If  03o  is  the  same  starting  speed  and  o>'  the  new 

speed  at  time  t, 

,  -Rt/o" 

03  =  03*  e  (£>) 

From  equation  (5>  and  equation  C6>  one  may  deduce  that 

1,  _ Zll-S: _  or  that  J  -  AJ  ^  ^  (7) 

’  U  &  U,  ml 

o>  Ui 

where:  01* is  the  speed  of  the  rotating  body  at  the  moment  all  supplies 
of  power  are  removed. 

o>'  is  the  speed  at  time  t  with  added  inertia, 

o)  is  the  speed  at  time  t  without  added  inertia. 


at  is  the  added  moment  of  inertia. 
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Similarly  one  may  deduce  that 


R  _  ' 

i  (o 


AJ 


(Jo 
CO  ' 


( 8 > 


Equation  (4)  can  be  modified  to  allow  for  a  purely  frictional  torque 
as  well  as  damping  torques.  If  Tf  is  the  constant  frictional  torque, 

the  new  differential  equation  is 

J  to  +  R  co  +  Tf  .  o  (9) 


With  the  use  of  Laplace  Transforms 

J" 6l>  ■+  Rca  4  Tf  _  O  oecames 
J  {  S  CJ  ( i)  -  }  4  R  it)  C  s)  +Tf_ 


(9) 


Thus 
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s  (s.f^ 


(io> 


Using  Inverse  Laplace  Transforms 


-Rt/j  -Rt  /cr 

(A  U)  -  CA0  e  -jr  C  i  -  e  ) 
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If  the  moment  of  inertia  can  be  increased  without  changing  R  or  TF 
then 

Rt/J’ 


_  ]V  (  i 
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(A  Lt)  -  (A0  e 

From  equation  (11) 

alt)  +  I f  -  ( 6)0  +  TF  )  e 

From  equation  (12) 
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Rearranging  equation  (13) 
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or 


R 


J  { 

^  uCt)  +  _Z>  J 
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Similarly  equation  (14)  gives 
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Since  T  -  J~  +  AJ~  equation  (17)  is  equivalent  to 

v~  j 

0~  _  A  J  (  o 


Y«0  4  Tr/f,  j 


4)' (it)  +  Tr/R 


l  W(f)  +  Tf/B 

Substituting  into  equation  ( 13A) 
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R  = 


J_  aJ 
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f  -  *  F/R  -  [ 

(  to1  (t)  4  T>  /R  J 

f  6)70  4  Tp/R  | 

l  toU)  +  TF  /R  J 
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+  ip/r 


O)  (fc )  +  tf  /p  . 


(17) 
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where  6)„  is  the  speed  of  the  rotating  mass  when  all  power  sources 


are  disconnected 

6)(t)is  the  speed  of  the  rotating  mass  when  it  is  coupled  to 
another  rotating  mass  and  the  time  is  t  seconds  after  all 
power  sources  are  disconnected 
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cjc-t)  is  the  speed  of  the  rotating  mass,  alone,  t  seconds  after 
all  power  sources  are  disconnected. 

At  is  the  added  moment  of  inertia 
R  is  the  constant  damping  function. 


Ordinarily  Friction  forces  and  torques  are  unknown  and  R 
is  distinctly  unknown.  Consequently  the  term  tf/r  is  unknown. 


If  a  different  AJ,  say  AJ'  ,  is  employed  and  the  reasoning 
repeated  then  one  must  obtain  the  same  value  of  J  as  in  equation 
(18).  Thus 


J  = 


aj'. 
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C0o  +  Tp  /  R 
6>"<t)  +  TF/~K 

u>"U)  *  Tr/B 
(A  (t)  -t  Tp  /R 
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where  Ajis  the  new  increase  in  moment  of  inertia 

(od)is  the  new  speed  at  time  t  with  the  new  increase  in  moment 
of  inertia 

colt) is  the  speed  of  the  rotating  mass  at  time  t  without  any 
any  increase  in  moment  of  inertia. 
ca0  is  the  speed  of  the  rotating  mass  the  instant  all  power 
sources  are  removed. 
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Clearly  an  iterative  procedure  could  solve  equations  ( 18A) 
and  ( 19A)  for7>/f?and  R  could  be  calculated  from  either  equation. 

Summary : 

(1)  The  standard  technique  for  determining  the  polar 
moment  of  inertia  tends  to  indicate  a  polar  moment  of  inertia  which 
is  less  than  the  actual  value  because  it  neglects  the  purely 
frictional  torques. 

(2)  Modifications  which  allow  for  pure  friction  torque 
result  in  formulae  which  require  iterative  procedures  for  their 
solution. 


The  writer's  technique  for  determining  the  polar  moment 
of  inertia  does  not  give  the  accuracy  of  the  retardation  technique 
but  it  does  have  the  advantage  that  it  may  be  built  in-situo  from 
scrap  material  and  for  most  practical  purposes  the  accuracy  is 
sufficient.  This  technique  is  outlined  in  Section  IV  of  this 
chapter. 


Since  the  polar  moment  of  inertia  and  damping  figure  so 
prominently  in  the  mechanical  equations  of  motion  a  detailed  study 
of  friction,  windage,  viscosity  and  the  polar  moment  of  inertia  was 


undertaken. 
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IV.  FRICTION  BY  MODIFIED  SUPPLY  OF  LOSSES  TECHNIQUE 

The  wound  rotor  machine  studied  was  normally  coupled  to  a 
compound-wound  direct  current  machine,  the  armature  and  all  field 
connections  of  which  were  separately  accessible.  The  shunt  field 
of  the  direct  current  machine  was  separately  excited  from  a  constant 
current  supply  while  the  series  fields  were  unused. 

Combined  armature  brush  and  copper  less  was  determined  at 
ambient  temperature.  Since  armature  currents  during  the  test  were 
low  no  temperature  corrections  were  applied. 

For  several  values  of  decreasing  speed,  speed,  current 
and  armature  terminal  voltage  were  recorded  as  in  the  conventional 
"Efficiency  by  Supply  of  Losses"  Method,  this  being  done  for  the  case 
when  the  wound  rotor  was  coupled  to  the  direct  current  machine.  When 
the  rotor  stalled  the  wound  rotor  induction  machine  was  uncoupled  and 
the  mechanical  losses  of  the  direct  current  machine  were  used  as  the 
load  for  the  direct  current  machine. 

The  gross  power,  PT  ,  applied  to  the  terminals  of  the 
armature  of  the  direct  current  machine  is  composed  of  the  following 
terms : 

(a)  ,  the  hysteresis  power  loss  in  the  armature,  which 

/■e 

is  proportional  to  <p  f)  where  <p  is  the  flux  in  the  armature  and  ^ 
is  the  velocity  of  the  armature. 

(b)  PE  ,  the  eddy  current  power  loss  in  the  armature, 
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2 

which  is  proportional  to  n  . 

(c)  PR  ,  the  resistance  power  loss  in  the  armature,  which 
depends  on  the  brush  and  copper  resistance  and  the  square  of  the 
armature  current. 

(d)  PF+W  ,  the  total  friction  and  windage  power  loss  of  the 

armature . 


If  one  could  keep  the  armature  current  sufficiently  small 
then  the  armature  flux  would  be  independent  of  armature  current.  If 
one  could  assume  that  radiation  is  the  chief  manner  in  which  the 
armature  is  cooled  then  the  value  of  PR  would  be  independent  of 
speed.  If  one  could  assume  that  friction  and  windage  dominate 
viscosity  forces  then  one  could  and  would  assume  that  '"V+w  is 
independent  of  temperature.  Unfortunately  none  of  these  assumptions 
is  valid  but  all  of  these  assumptions  are  necessary  for  the  analysis 
to  follow. 


If  the  subscript  A  is  used  to  denote  a  study  when  the  wound 
rotor  is  coupled  to  the  direct  current  machine,  and  subscript  B  is 
used  to  denote  that  only  the  direct  current  machine  is  running,  then 


pha  +  Pea  +  Pra  +  pf,wa  =  Pta 


or 


PHa  +  PEa  +  PF  t  WA  =  PTa 
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Subtracting  equation  (2)  from  equation  ( 1) 


<P«  *  PE.)  '  (  P»  '  P0  +  (FW  '  P^~»)  =  <  PrA  -  PR„)  -  <  PT  - 


(3) 


Now  all  terms  on  the  left  hand  side  of  equation  (3)  are  flux  and/or 
speed  dependent.  If  the  independent  variable  is  speed  and  the  flux 
is  constant  certain  terms  in  equation  (3)  cancel  and  thus  the  friction 
and  windage  power  in  the  alternating  current  machine  must  be  given  by 
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(4) 


The  exact  procedure  involved  will  be  fairly  obvious  from  a  study  of 
the  several  figures  in  Appendix  C. 


V.  MECHANICAL  MEASUREMENT  OF  MECHANICAL  DAMPING 

The  results  obtained  by  the  technique  of  Section  IV  were 
not  completely  reliable  because  there  seemed  to  be  no  way  of 
reproducing  them  at  later  dates.  Consequently  this  writer  felt 
that  an  attempt  should  be  made  to  determine,  mechanically,  the 
mechanical  damping  and  polar  moment  of  inertia. 

This  technique  involved  wrapping  aluminum  wire  around 
a  shaft,  passing  it  over  a  pulley  and  applying  a  large  weight  on 
the  other  end  of  the  wire. 

The  technique  necessitated  the  assumptions  as  follows: 

(a)  The  pulley  had  zero  moment  of  inertia 

(b)  The  wire  was  weightless 
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(c)  The  wire  was  inextensible. 


If  W  is  the  applied  force  and  r  is  the  shaft  mean  radius, 
then  wr  is  the  applied  torque.  However  a  purely  frictional  torque 
can  be  deduced  from  fig.  4  Appendix  C  and  deducted  from  Wr  . 

Unfortunately  the  solution  of  the  differential  equation, 
which  describes  the  motion  of  the  weight,  is  cumbersome  and  involves 
the  solution  of  a  transcendental  equation. 


The  application  of  Newton's  Third  Law  gives 
T  60  +  R  a)  =  w'r 

Applying  Laplace  Transforms  to  equation  (1)  gives 
T(sl0s  -  50o-  0.)  +  R  (50J-0.)  .  wr 

Therefore 
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If  t,  represents  the  time  when  the  weight  reaches  minimum  altitude  the 


Inverse  Laplace  Transform  of  (2)  gives 
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If  the  weight  starts,  at  height  h,  ,  from  rest  and  proceeds 
downwards,  -  h.  where  h  is  the  distance  the  weight  drops  and  r  is 


the  mean  radius.  Thus 
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After  the  weight  reaches  its  minimum  altitude  the  applied 
torque  on  the  shaft  reverses.  The  velocity  at  minimum  altitude  is 
given  by  equation  (4)  and  the  value  of  ree  is  given  by  equation  (3). 
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When  the  weight  reaches  its  maximum  height, 
thus 
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Resubstituting  for  (  t-  e  ) in  equation  (5)  reveals  that: 
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While  equation  (9)  appears  to  be  very  interesting  it  is  of  little  use 
since  ts  ,  and  tt  ,  are  unknown  quantities.  Furthermore,  there  does 
not  seem  to  be  any  method  to  determine  -e,  ,  and  t‘  ,  even  though 

t  *  t,  +  t]  can  be  measured  and  h  and  h'  can  be  measured.  The  only 
truly  satisfactory  procedure  would  be  to  measure  t,  and  t'  separately. 
This  could  be  done  with  two  photo-electric  cells.  If  so  done  the 
damping  function  R  could  be  determined  to  very  great  precision.  Then 
substitution  of  R  into  equation  (3)  would  determine  J  in  terms  of  a 
transcendental  equation. 


Clearly  an  iterative  procedure  is  involved  in  the  above 
solution  but  if  the  same  care  is  employed  as  in  the  retardation 
technique,  excellent  accuracy  would  be  expected. 

In  order  tc  achieve  a  solution  to  the  problem,  the 

rr  • 

differential  equation  was  assumed  to  be  Jco  +  R  a>  =  -  J  wr  to  dt 
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which  has  the  solution 

_  zt 

9  =  A  e  23~  cas  t*/(  &T-/JL\X) 

V  J  \2j/ 

the  period  of  which  is 


r  = 


l  7T 


The  nature  of  the  errors  in  the  above  approximation  is  unknown  and 


cannot  be  predicted. 


A  very  quick  method  for  determining  the  moment  of  inertia 
has  been  demonstrated  and  formulae  have  been  derived  to  indicate 
possible  improvements. 


Appendix  C  explains  the  procedure  employed  to  obtain  the 
polar  moment  of  inertia. 


VI.  MOMENT  OF  INERTIA  BY  ELECTRICA1  MEANS 


If  the  direct  current  machine  with  constant  field  current 
is  used  to  accelerate  the  rotor  of  the  alternating  current  machine 
then  theory  indicates  as  follows: 


1  J  (  ez  -  ©  *  )  -  f  vi  dt  _  [CPHdt_  '  Pe  dt  _  rSv+L0 dt  _  f£p„  dt 

2  v  ^  y  ..  t.  q  y  ^  W 


where 


0  is  the  angular  velocity  at  time  t 
60  is  the  angular  velocity  at  time  t  =  0 


V  is  the  armature  terminal  voltage 


P  p  P  p 

rH  »  £>  >  rR 


are  as  defined  in  Section  IV. 


Clearly  PF+W 


and  PR  are  reasonably  well  known.  However, 


. 


: 
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P*  and  PE  are  not  known  and  one  requires  a  "Separation  of  Losses" 
method  to  determine  them. 

In  order  to  ascertain  whether  or  not  the  idea  might  be 
fruitful,  several  attempts  were  made  to  correlate  armature  terminal 
voltage  with  speed.  No  correlation  could  be  found  and  the  attempts 
to  measure  PH  and  PE  were  not  made. 

Several  factors  explain  the  failure  to  obtain  useful 
results  from  this  method: 

(a)  Armature  reaction  reduces  the  generated  back  voltage. 

(b)  Hysteresis  and  eddy  current  power  are  dependent  upon 
flux  density  and  speed  and  the  generated  back  voltage  is,  therefore, 
partly  dependent  on  these  quantities. 

(c)  Self-inductances  in  the  shunt  field  prevent  the 
reestablishment  of  shunt  field  flux. 

This  method  could  be  improved  by  using  a  cathode  ray 
oscilloscope  equipped  with  an  electronic  switch  to  display 
instantaneous  speed  from  a  tachometer  generator. 
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CHAPTER  V 


METHOD  OF  SOLUTION  AND  CONCLUSIONS 

I.  METHOD 

Rough  estimates  indicate  that  a  time  increment  of  one  milli¬ 
second  would  be  too  large  and.  thus  a  large  volume  of  computer  output 
would  be  expected.  The  choice  of  digital  computer,  LGP  30,  was  un¬ 
fortunate  in  several  respects:  firstly  the  outputting  devices  for 
this  computer  are  very  slow;  secondly  fixed  point  programming  was 
necessary  in  order  to  gain  internal  speed  of  the  computer;  thirdly 
fixed  point  programming  restricts  the  range  of  variables  and  finally 
errors  are  more  frequently  inserted  in  fixed  point  programming.  Hie 
other  digital  computer,  IBM  1620,  had  a  heavy  work  load  and  was  virtually 
unavailable  to  this  writer.  Since  fixed  point  programming  was  necess¬ 
ary,  a  detailed  knowledge  of  the  parameters  was  necessary  before 
programming  could  be  attempted. 

The  transient  currents,  rotor  accelerations,  rotor 
velocities  and  rotor  displacements  are  considered  from  the  moment 
two  systems,  at  slightly  different  frequencies,  are  paralleled  through  the 
machine.  The  same  phenomena  are  studied  when  both  systems  experience  a 
line-to-line  fault,  at  the  machine,  at  the  moment  of  paralleling.  In 
both  cases  it  was  assumed  that  stator  currents  had  achieved  steady  state 
before  paralleling  the  rotor  circuit.  Stator  and  rotor  frequencies  were 
376.991119  and  350.000  electrical  radians  per  second,  respectively, 
initial  rotor  velocity  and  displacement  were  0.00  electrical  radians 
per  second  and  0.0192  electrical  radians  respectively. 
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The  writer  attempted  to  program  the  solution  directly  from 
the  differential  equations  found  in  Chapters  I,  II  and  III.  The 
results  of  such  a  program  were  most  disappointing .  These  equations 
were  then  rearranged  as  in  Table  I  of  this  chapter  and  a  new  program 
was  written.  This  program  produced  suitable  results. 

The  time  increment  used  was  0.2  milliseconds.  This  time 
increment  appears  to  be  almost  optimal  and  should  not  be  changed. 

Fig.  1  displays  the  typical  counter-clockwise  motion  with 
increasing  time  and  therefore  the  locus  is  that  of  a  positive  sequence 
function.  The  locus  is  almost  hypocyc loidal  in  nature  and  displays 
marked  periodicity  in  both  phase  and  magnitude.  The  peak  magnitude 
is  approximately  6.5  amperes.  The  locus  is  of  the  positive  sequence 
stator  current. 

Fig.  2  displays  the  typical  counter-clockwise  motion  with 
increasing  time  and  therefore  the  locus  is  that  of  a  positive  sequence 
function.  As  in  Fig.  1  the  locus  is  almost  hypocyc loidal  and  displays 
marked  periodicity  in  both  phase  and  magnitude.  The  peak  magnitude 
is  approximately  13.5  amperes.  The  locus  is  of  the  positive  sequence 
rotor  current. 

Figures  3  and  4  display  the  rotor  angular  accelerations  in 
electrical  radians  per  second  per  second.  The  accelerations  are 
periodic  in  nature  passing  through  the  time  axis  at  roughly  0.0166, 
0.0330,  0.050,  0.0666,  0.0833,  0.0995,  0.1160  etc.  seconds.  These 
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times  coincide  in  a  rough  manner  with  the  periods  of  the  stator 
voltages.  The  peak  acceleration  observed  was  about  370  electrical 
radians  per  second  per  second. 

Figures  3,  6  and  7  display  the  rotor  angular  velocity  in 
electrical  radians  per  second.  The  velocity  is  a  complex,  periodic, 
and  ill  behaved  function  which  is  clearly  unstable. 

Figures  8,  9  and  10  display  the  rotor  angular  displacement 
in  electrical  radians.  The  displacement  function  clearly  displays 
instability . 

The  digital  computer  was  programmed  in  such  a  manner  that 
a  fault  could  be  imposed  at  any  time  from  start  of  double  excitation 
until  fifteen  seconds  after  the  start  of  double  excitation.  The 
following  types  of  faults  are  considered: 

(a)  A  three  phase  fault  on  either  stator  or  rotor 

(b)  A  line- to-ground  fault  on  either  stator  or  rotor 

(c)  A  line-to-line  fault  on  either  stator  or  rotor 

(d)  A  line-to-line  fault  on  both  stator  and  rotor. 

The  three  phase  fault  is  very  uninteresting  insofar  as  the 
doubly  excited  induction  machine  is  concerned.  A  severe  but  brief 
power  inrush  occurs  to  cause  the  rotor  to  accelerate  to  normal  speed. 
The  faulted  system  is  clearly  in  serious  trouble  but  the  unfauited 
system  and  the  doubly  excited  induction  machine  are  not  seriously 


affected . 
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The  line- to-ground  fault  could  pass  zero  sequence  current 
through  the  machine  but,  with  the  neutral  points  ungrounded,  zero 
sequence  currents  would  be  negligible. 

The  line- to- line  fault  on  either  stator  or  rotor  applies 
negative  sequence  potentials  to  the  side  faulted,  causing  negative 
sequence  currents  to  flow  in  the  side  faulted,  inducing  potentials 
of  negative  sequence  in  the  unfaulted  side  and  causing  currents  of 
negative  sequence  to  flow  in  the  unfaulted  side.  Clearly,  however, 
the  negative  sequence  current  in  the  unfaulted  side  would  have  some 
impedance  through  which  to  flow  and  this  voltage  drop  would  be 
transformed  to  the  faulted  side.  Thus  the  fault  current  in  the 
doubly  excited  induction  machine  would  be  limited  to  some  extent. 

The  line- to- line  fault  on  either  stator  or  rotor  suggests 
an  examination  of  the  case  where  both  stator  and  rotor  have  a 
line-to-line  fault.  The  negative  sequence  currents  would  be  most 
severe  under  these  circumstances. 

Clearly,  if  the  doubly-excited  wound  rotor  induction  machine 
can  withstand  the  worst  of  four  cases  it  can  withstand  all  cases. 
Clearly  the  wisest  choice  is  to  study  the  worst  possible  case  as  far 
as  fault  currents  are  concerned.  For  this  reason  the  program  was 
designed  so  that  symmetrical  double  excitation  or  line-to-line  faults 
on  both  rotor  and  stator  could  be  selected  at  the  discretion  of  the 


operator. 
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Figure  11  displays  typical  counter-clockwise  motion  and  is 
therefore  the  locus  of  a  positive  sequence  function.  It  is  the  locus 
of  the  positive  sequence  rotor  current  under  the  fault  condition 
selected  above.  Fig.  11  is  quite  similar  in  nature  to  Fig.  2. 

Fig.  12  is  also  the  locus  of  a  positive  sequence  function. 

It  is  the  locus  of  the  positive  sequence  stator  current  under  the 
fault  conditions  described  above.  Fig.  12  is  very  similar  to  Fig.  1. 

Fig.  13  displays  the  negative  sequence  stator  current  under 
the  fault  conditions  described.  The  peak  negative  sequence  current 
is  about  9.8  amperes.  Irregularities  are  due  to  a  slight  error  in 
the  program  which  had  a  tendency  to  neglect  the  imaginary  component 
of  rotor  negative  sequence  resistive  potential  drop.  This  error  was 
discovered  after  all  the  loci  had  been  plotted.  Changes  were  made  to 
the  program  but  the  program  was  not  re-run. 

Fig.  14  displays  the  negative  sequence  rotor  current.  The 
peak  magnitude  of  negative  sequence  rotor  current  is  about  13  amperes. 
Any  irregularities  are  due  to  the  error  described  under  Fig.  13. 

Fig.  15  displays  the  rotor  acceleration  in  electrical  radians 
per  second  per  second.  The  peak  acceleration  is  about  -  520  electrical 
radians  per  second  per  second.  The  rotor  acceleration  is  periodic  and 
has  nearly  the  same  periods  as  in  the  unfaulted  case.  However  the 
accelerations  are  quite  different  in  the  two  cases. 


Fig.  16  displays  the  rotor  velocity  in  electrical  radians 
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per  second  under  the  fault  conditions.  The  velocity  is  in  the  wrong 
direction  and  therefore  indicates  instability. 

Fig.  17  displays  the  rotor  displacement  in  electrical  radians 
under  the  fault  conditions.  The  displacement  is  in  the  wrong  direction 
and  indicates  instability. 

Figures  1  through  17  were  the  lesults  of  calculations 
performed  by  a  digital  computer. 

A  test  of  the  doubly-excited  induction  "motor"  was  made. 

With  a  small  shaft  load  the  writer  found  that  two  synchronous  systems 
could  be  asynchronously  connected  provided  that  they  were 
"asynchronously"  connected  at  the  same  frequency  before  any  changes 
in  frequency  occurred.  Even  so,  large  fluctuations  in  currents  and 
powers  nullified  the  test  except  to  indicate  as  follows: 

(a)  Power  can  be  asynchronously  interchanged 

(b)  The  magnitude  of  power  flow  can  be  changed  by  changing 
the  ratios  of  stator  to  rotor  potentials. 

When  the  small  shaft  load  was  uncoupled  the  doubly  excited 
induction  machine  became  unstable  and  the  experiment  w7as  stopped. 
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II.  CONCLUSIONS 

For  the  machine  studied  in  the  laboratory,  the  digital 
computer  study  indicates  instability  if  the  difference  in  applied 
frequency  exceeds  26.991118  radians  per  second.  This  is  equivalent 
to  some  4.3  cycles  per  second.  Probably  instability  would  occur 
at  lower  speeds  and  frequency  differences  although  such  tests  and 
calculations  were  not  made. 

This  instability  is  attributed  to  lack  of  rotor  damping 
circuits.  A  digital  computer  study  of  such  machinery  with  rotor 
damping  circuits  would  require  a  digital  computer  of  much  greater 
capacity  than  the  Royal  McBee  LGF  30. 
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22 

23 
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x  U  /  3  35  4- 
X  tA  3?  C  O  3 
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STOP 
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I- 
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I 
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1 
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r 
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1 
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r 
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i 
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r 
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f 

1 -  1 
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^9 

1 
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1 
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■ 
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. 
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— — — j — — - 
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X 

■ 
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i 

1  1  1  1  1  1  1 
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i 

i  i  i  |  i  i  r 

i  ,  46 

/  5  S  O  6  / 

, 

1  1  1  1  1  1  1 
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’ 
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1  1  1  1  1  1  1 

■  .  4,8 
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1  1  1  1  1  1  1 
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f 
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i 

• 
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t 
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i 

X 

L 
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i 
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» 

@3 
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1  54 

1  l^|5|/|/|0,5 

i 
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55 
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L  |  i  i  |  i  i  i.1 

X 

_ 
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» 
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r — 1 — 1 — 1 — 1 — L 
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i 
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i 
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L 

1  1  1  1  1  1  1 
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i  i  i  i  i  l  1 

i 
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i 

- L - 1_ 1_ 1_ 1_ 1_ 1 - 
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6J 
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i 
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’ 
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C941B3  “CSIE  TORONTO  C»N»0» 

CONDITIONAL  STOP  CODE 


CARRIAGE  RETURN 


The  MeBee  Company  Limited 


JAN  57 


LG  P-30  CODING  SHEET 

NO.  1 


PAGE  .  // 


40 


OF 


PROG  NO. 


PREP  BY 


CKD  BY 


DATE 


TRACK 


Lgp-jO  cob! 

PAGE  lZ  _  OF 

NO  2 


PROG  NO. 


PREP  BY 


CK'D.  BY 


DATE 


M 


TRACK 


PROGRAM  INPUT 

CODES 

STOP 

LOCATION 

INSTRUCTION 

OP.  ADDRESS 

STOP 

CONTENTS  OF 

ADDRESS 

NOTES 

™  l 

1  1  -  J  1_ 1_ 1_ 1 

1 
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• 

—  i 

1  _l _ J _ 1 _ 1 - 1 - L 

1 

1 1  1  1  1  1  1— 

1  +j 3|2 

«  T 

_ l  i  a  i  U 4  /  j  A  /  3 

i 

| 

i  J  1  L  i  1_ L  . 

l_L_L_M 

L _ i _ i _ | _ i _ i  i  i 

,  /;  C|  /  |  O  & 

i 

<p  c ?/2 

/.y  • 

«  1 

1  1  -J | 1_ 1 ! 

i  ■  3,4 

l 

I 

1  1  J  1_ 1_ 1_ 1 

i  i  3,5 

i  x  i  U  |  /  4-  5  2 

X 

j 

1111  1  1  1 

■  1  1  3;6 

i  A  .  4^1 1  ^  j  3  i  /  i  / 

i 

% 

1 

1  1  1-1_ L l J_ 1 

i  1  3,7 

i  i  x  \C  1 2  |3  i  !  \£> 

l 

* 

. 

r  i 

1  1  1  |  1  1_ L 

,  ,  3,8 

i  /  , 3 , ^|3 

l 

diOf 

Ax 

1! 1 |  1  1  1 

i  L-3,9 

.X.y|/,4|4  2 
,  ,  /  ,A  ,  /  ,Z,S,S 

X 

1  1  1  |  1  1 ! 

1 

r&s/  psrf  <1  -  3Z  1 

k/  e  -*>  j  (S/s 

I  1  1  1  |  1  1  1 

,  ,  41 

1 

i  i  U  j_/  aA'  .4  & 

• 

1  1  1  |  1  1  1 

i  i  4j2_ 

,  X  /tt.  5,0,  O  7 

t 

real  part  & 

9 

,  i  4,3 

_jl  i  /  0, 6  3  <o  Z 

’ 

X 

.  1  1  1  |  1  1  1 

j _ l  4  4 

,  ,4,  /  7 

3  O  0|C>,  0,0  c 

1 

1  1  1  1  1  1  1 

i  ,  4,5 

l 

/  CB4- 

r —  — \ 

i  i  46 

2  0,0  O  0,0,  <o  <o 

l 

• 

-i  i  i  1  i  i  i 

,  ,  47 

1  1  1  1  1  1  1 

I 

a. 

i  i  i  i  i  i  i 

_l  i  4.8 

1  /  ,3.3.4 

X  A  y! ±3 ,0  v7 

i 

/  @4-  ■ 

,  4l9 

I 

l 

real  parX  J  /- 

■^za./e'y^  ( 5/3 

i  i  1  1  i  i  i 

i  i  i  i  i  i  i 

j_  i  ^0 

X  M ,  /  ,3  0  & 

— 1 — 1 — ! — 1 — J — L 

,  5,1 

r 
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l 

X  _ 

i  i*i  |  i  i  i 

,  i  5,2 

— —  - 

X  S.SOJ  O 
i  j  X  U  |  /  (o  !  4- 

t 

i  i  i  |  i  i  i 

i  i  i  t  i  i  i 

,  i  5,3 

i  i  i  i  i  i  i 

,V  i  5,4 

,  x  XxN\\S\Cf\t  kZ_ 
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_ I _ i»X  rU  |  1  7  3  2 
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1 

* _ 
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I 

• 
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l 

i 

- 1 - 
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— * — 1 — -i — -* — — 1  — 

,  5,9 

X- 

✓ 
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r 
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■  — 
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*  1 - 1 - L_  -L  J _ l  -  L 

61 

x  U  \l  i-3~ i  / 

— — * - 1 - 1_ L  J_ 1_ 1 - 

1 _ L_ - 1 - 

.i  ,.l  6  2 

— i — i — ir-  |.'.-i — 
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• 
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J  L 
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-h 
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33 


35 
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,  ,  X  U  |  l  S  3  4 

x  /vi ,  /  ,  2  S  6 
A  ll  \  5  3  8 


i 


36 


J L 


37 


4- 1  Q\  Oj  O i  0\0\0 

L .  1 


J _ L 


38 


x  A  l  2  6  O 


J  L 


39 


X  /Vf  /  2  2  5 


40 


L L 


!  5  4  / 


i  i  **i 


41 


i  4,2 
,  .  43 

.  ■  4,4 


X  5  <2>  0  6  3 


j _ l  X  U  yl  \5 3 

_l  i  X  jM\  LlL  i  QaL 


_J _ L 


45 


j 


_ i _ 


i  l  48 


i  i  i  4,9 

»  i  5,0 

■  ■  5,1_ 

i  i  5j2_ 


L 


,5,3 
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61 
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63 
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i  i  X  iCj  /  i<3  /  3 
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- 


m 
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K 
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03 
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O 

.  0  8 
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OP 
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J _ L 
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J 

1_ 


+- 
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-1 — i - i - L 


J _ L 
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-1  L- f 


_1. 
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_L 
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J _ L 
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4- 

-4- 
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-h 
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1 

J _ L- 
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f 
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| 
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t 

} 
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AEPENDIX  B 


OPERATING  PROCEDURE 


Stage  A: 

Before  depressing  the  "Power  On"  button  the  following  checks 
should  be  made: 

(1)  The  "input  selector"  switch  should  be  to  "typewriter", 

(2)  The  "output  selector"  switch  should  be  to  "  typewriter'^ 

(3)  "Reader  Power"  should  be  up. 

(4)  "Punch  Power"  should  be  up. 

(5)  On  the  consul 

(a)  the  "manual"  button  should  be  down. 

(b)  the  button  between  "Standby"  and  "Operate" 
must  be  depressed. 

(c)  all  breakpoints,  the  six  bit  input  switch  and 
the  transfer  control  switch  must  be  up. 

(d)  the  "execute"  button  must  be  up.  The  LGP  30 
at  University  of  Alberta  has  a  tendency  for 
this  button  to  stick. 

Stage  B: 

After  the  checks  of  Stage  A  have  been  completed  the  operator 
may  depress  "power  on"  switch.  When  the  "operate"  light  appears  the 
operator  may  proceed  to  Stage  C. 

Stage  C: 

(1)  The  operator  should  turn  flexowriter  power  "on"  and 
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the  "reader  power"  button  should  be  depressed. 

(2)  The  flexowriter  connector  switch  should  be  off  and 
the  manual  input  lever  on  the  flexowriter  should  be  down.  All  other 
levers  on  the  flexowriter  should  be  up. 

(3)  Paper  tape  number  1  should  be  inserted  in  the 
flexowriter  mechanical  reader. 

(4)  The  operator  should  depress  the  "Start  Read"  lever. 

This  will  result  in  C3Q00 '  being  typed  and  being  read  into  the 
accumulator. 

(5)  The  operator  should  depress  "Fill  Instruction"  on 
consul.  This  would  result  in  an  order  being  placed  into  the 
instruction  register  namely  an  order  to  clear  the  contents  of  the 
accumulator  into  track  62  sector  00. 

(6)  The  operator  should  depress  "start  read".  The  result 
would  be  I00J8  being  typed  and  read  into  the  accumulator. 

(7)  The  operator  should  depress  "one  operation",  "execute 
instruction",  and  "manual"  buttons  in  that  order.  The  "one  operation" 
puts  the  computer  into  a  mode  where  the  "execute  instruction"  can  be 
effective  and  the  instruction  filled  in  stage  C(4)  clears  an 
instruction  (input)  into  track  62  sector  00. 

(8)  Stages  C(4),  C(5),  C(6)  and  C(7)  should  be  repeated 
until  the  instruction  U3Q00 '  has  been  loaded.  The  instruction  which 
loads  U3Q00 '  is  H3Q10'  which  means  store  the  contents  of  the  cumulator 
( U3Q00 ' )  in  track  62  sector  04.  The  accumulator  will  retain  U3Q00. 

(9)  Depressing  "Fill  Instruction"  will  place  the  instruction 
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U3Q00  into  the  instruction  register  and  such  will  be  necessary. 

(10)  The  operator  should  depress  "one  operation",  "execute 
instruction"  and  "manual". 

(11)  The  operator  should  depress  "Reader  Stop"  and  insert 
tape  number  2  into  the  photoreader,  turn  the  "Input  Selector  Switch" 
to  "photoreader". 

(12)  The  operator  should  depress  "One  Operation",  "Normal" 
and  "Start"  on  the  consul.  The  result  should  be  that  the  computer 
should  load  a  paper  tape  program. 

(13)  The  computer  will  stop  without  having  taken  all  of  the 
program.  The  remainder  will  be  composed  of  two  parts;  a  portion,  in 
hexadecimal,  self  loading  in  track  06  and  a  portion  in  binary-coded 
decimal  which  will  be  relocatable. 

(14)  When  the  operator  presses  "start"  on  the  consul  the 
computer  will  load  more  program  into  track  06  and  halt.  The  remainder 
of  the  tape  should  be  of  no  consequence  to  this  paper. 

(15)  The  operator  should  insert  tape  number  three  into  the 
photoreader  and  depress  the  "start"  button.  The  computer  will  stop 
near  the  end  of  the  tape.  If  the  operator  should  desire  to  study  the 
doubly-excited  induction  machine  with  artificial  damping  torques  he 
should  remove  the  program  tape.  If  he  wishes  to  study  the  doubly- 
excited  induction  machine  without  artificial  he  should  depress  "start". 

(16)  The  operator  should  insert  an  appropriate  data  tape  in 
the  photoreader  and  press  "start". 

(17)  The  operator  should  insert  a  quantity  of  paper  tape  in 
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the  punch,  depress  all  break  points  and,  if  a  symmetrical  excitation 
is  to  be  studied,  the  transfer  control  button  should  also  be  depressed. 

(18)  The  last  "word"  on  the  data  tape  should  be  .0000942'. 

The  punch  power  should  be  on.  The  operator  should  then  depress  the 
"start11  button.  After  the  computer  commences  punching,  the  operator 
should  turn  the  input  selector  switch  to  typewriter.  Then  the  operator 
should  turn  off  the  "reader  power"  and  the  flexowriter. 

Stage  D: 

(1)  If  the  computer  must  be  released  to  others  the  operator 
should  wait  until  the  computer  starts  punching  and  then  lift  bree.k 
point  8.  He  should  depress  "Manual  Input"  lever  on  flexowriter. 

(2)  When  the  computer  halts,  the  operator  should  turn  on 
the  flexowriter  power,  insert  tape  number  four  into  the  mechanical 
reader  and  depress  "One  Operation",  "Clear  Counter",  "Normal",  "Start". 
When  the  light  appears  on  the  f lexowriter ,  the  operator  should  depress 
"tape  feed"  and  tear  off  the  data  output.  He  should  depress  "tape  feed" 
again  and  lift  the  manual  input  lever.  This  will  result  in  a  retention, 
on  tape,  of  valuable  hexadecimal  data.  It  should  be  carefully  labelled. 

Stage  E: 

(1)  If  the  operator  wishes  to  continue,  at  a  later  date,  from 
DC  2 )  he  should  load  the  program  (tape  number  three)  by  depressing  "one 
operation",  "clear  counter",  "normal,  "start"  after  tape  number  three 
has  been  inserted  in  the  photoreader  and  the  "input  selector  switch"  has 
been  turned  to  "reader".  The  same  procedure  should  be  used  to  load  the 
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hexadecimal  data  of  D(2). 

(2)  The  operator  should  turn  the  punch  power  on  and  turn 
the  "input  selector"  switch  to  "typewriter".  When  the  flexowriter 
light  appears  after  "one  operation",  "clear  counter",  "normal",  "start", 
the  operator  should  type  .0002445' .  He  should  depress  all  break  points 
and  the  "transfer  control"  button  should  be  as  in  D(2).  Depressing 
"Start"  button  will  cause  the  computer  to  continue  from  D(2). 
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SAMPLE  DATA  INPUT 


; 0000951 ’+0110003 ’072900 ’000’ * 

; 0001029 ’ +0150014 ’018820 ’ 000 ’ - ’ 

; 0001031 ’ +0110003 ’ 591000 ’ 000  ” 

,* 0001039 ’ +0110003 ’ 591000 ’ 000 ” 

; 0001100 ’ +0130014 ’ 081490 ’ 000  ” 

; 0001130 ’ +0110003 ’ 473000 ’ 000 ’ - ’ 

,* 0001132 ’ +0110003 ’ 473000 ’ 000  ” 

,* 0001153 ’ +0110001 ’ 021100 ’ 000 ’ - ’ 

,* 0001201 ’ +0110001 ’ 137600 ’ 000 ’ - ’ 

; 0001263 ’ +0110001 ’ 021100 ’ 000 ’ - ’ 

,* 0001315 ’ +0110003  *  591000 ’ 000  ” 

,* 000133 5  T +0110003 ’ 562000 ’ 000  ” 

,* 0001353 ’ +0110003 ’ 562000 ’000’ - ’ 

; OOOI363 ’ +0130014 ’ 018820 ’ 000 ’ - ’ 

; 0001505 ’+0110001’ 137600 ’000’ - ’ 
,*0001555’  +0110001  ’ 021100  ’  000  ’  -  ’ 

,*  000  i64o  ’  +0110003  ’  072900  ’  000  ’  -  ’ 

; 0001661 ’+0110001 ’ 137600’ 000 ’ - ’ 

; 0001700 ’ +0110003 ’ 562000 ’ 000  ” 

; 0001740 ’ +0110001 ’ 137600 ’ 000 ’ - ’ 

; 0001762 ’+0100000 ’ 986000 ’ 000 ’ ’ 

,*  00024 57 ’xz 0003  ’  ,*  0002 501 ’xz0003  ’ 


,*0005000’ +0700000  ’ 
383000 ’000” 
003371’000” 
049550’000” 

062000 ’000” 

000519 ’900” 
376991’ 119” 

350000 ’000’ ’ 
+0230009’ 

000123 ’000” 
006420’000*-’ 
,0000006’””” 
+0100000’ 

192000 ’000” 
+0430007’ 

000223 ’000” 

000149 ’000”””” 
+0230009’ 

128000 ’000” 
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MANUALLY  LOADED  BOOTSTRAP 
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Explanation  of  input  data. 


*,0000951 ' 
+0150014' 


xz0003' 


means  “start  filling  into  track  09  sector  51". 

means  "the  next  1  word,  9  digits  long  has  the  decimal  point 

after  the  fifth  (assumed)  character  and  is  to  be  multiplied 

_  i4 

by  2  before  storing".  Leading  zeros  need  not  be  punched. 
"Each  word  must  be  composed  of  three  parts:  the  most 
significant  6  digits,  the  next  3  significant  digits  and 
an  extra  stop  code  if  positive  or  -  if  negative", 
this  is  used  as  an  exponent.  All  exponents  are  stored  at 
q  =  29  in  the  subroutines  used  for  floating  point 
arithmetic.  xz0003  means  that  the  mantissa  as  a  decimal 
must  be  multiplied  by  1000.  A  negative  exponent  is  stored 
as  the  complement  of  the  positive  exponent  as  for  example 
;000243 1 ' ,000000 1 ' wwwwwwwj '  means  "Track  24  sector  31  is 
to  contain  the  hexadecimal  word  wwwwwwwj  which  is  -1  at 


q  =  29". 


Starting  in  track  50  sector  zero,  the  following  data  must  be  stored: 

r*,  ",  >,  Lp,  real  part  L«,  ,  imaginary 

part  ,  real  part  »  imaginary  part  iaz  , 

real  part  ,  imaginary  part  , 

initial  angular  displacement  from  stator  to  rotor  (radians) 

J,R. 

Starting  in  track  50  sector  20  the  stator  phase  voltage  must  be  stored. 
Starting  in  track  50  sector  22  the  rotor  phase  voltage  must  be  stored. 
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Starting  in  track  50  sector  59  the  difference  between  rotor  and  stator 

frequency  (radians/sec). 

Starting  in  track  24  sector  31  the  exponent  for  the  contents  of 

sector  32  track  50  sector  59. 


DATA  OUTPUT 

Three  rows  of  data  will  be  output: 

Row  #1  real  and  imaginary  parts,  positive  and  negative  sequence 

stator  currents. 

Row  #2  real  and  imaginary  parts,  positive  and  negative  sequence 

rotor  currents. 

Row  H 3  Angular  acceleration  of  rotor  electrical  radians  per 

second  per  second;  angular  velocity  of  rotor  electrical 
radians  per  second;  displacement  of  rotor  in  radians 
(a  floating  point  number)  as  for  example  .1919999810 
1  -  would  mean  .0191999981  radians;  electrical  displacement 
of  rotor  from  stator  in  radians  (a  floating  point  number); 
finally  the  time  elapsed  from  start  of  study,  in  seconds. 
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APPENDIX  C 


SUPPLEMENTARY  GRAPHS  AND  DATA 
I.  FRICTION  AND  WINDaGE  BY  ELECTRICAL  MEANS 

Figure  1  was  used  as  a  rapid  method  of  estimating  the  power 
lost  in  carbon  and  copper  when  the  shunt  wound  direct  current  machine 
was  running.  Table  I  displays  the  data  used  in  making  Fig.  4.  "Speed" 
was  measured  by  a  carefully  calibrated  Strobotac,  "potential"  was 
measured  across  the  armature  terminals  only,  "current"  was  measured 
by  a  miliivoltmeter  connected  across  a  carefully  matched  shunt. 

"Power"  was  immediately  calculated  as  the  product  of  current  and 
potential  and  was,  of  course,  the  gross  power  input  to  the  armature. 
"Copper  Power"  was  determined  from  Figure  1  and  subtracted  from 
"Power"  to  give  "Net  Power".  The  data  from  "Net  Power"  were  used  to 
produce  the  curves  marked  "A.C.  End  +  D.C.  End"  and  "D.C.  End  Alone" 
shown  in  Fig.  4.  The  difference  between  the  two  curves  was  used  to 
determine  the  curve  marked  "A.C.  End  Alone". 

Figures  2  and  3  were  attempts  to  gain  insight  as  to  the 
nature  of  friction  and  windage  of  bearings  and  to  exclude  brush 
friction  in  the  alternating  current  machine.  The  writer  could  only 
conclude  that  considerable  effort  would  be  needed  to  analyse  the 
phenomena  into  an  analytical  equation  and  that  such  an  equation 
would  not  be  of  much  use  in  any  case  since  the  results  were  not 
reproduceable  and  as  such  would  be  of  little  use. 
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The  following  phenomena  associated  with  friction  and  windage 
were  expected: 

(a)  a  portion  of  the  friction  and  windage  should  have  been 
linearly  speed  dependent. 

(b)  a  portion  of  the  friction  and  windage  should  have  been 
dependent  on  the  square  of  speed. 

From  equation  I  Chapter  III,  the  power  consumed  in 

,  2 

mechanical  damping  would  be  R (©)  but  this  quantity,  R,  would  be 

+2  -2 

measured  in  foct-lbs-sec  radians  whereas  the  quantities  that  could 

2  -2 

be  derived  from  Figs.  2,  3  and  4  would  be  watt-sec  revolutions  .  If 

the  quantities  could  be  derived  from  Fig.  4,  then  the  result  should  be 

2  -2 

multiplied  by  67.1  to  convert  watts-min  -rev  to  foot- lbs- sec-radians 

Figure  4  clearly  has  a  square  law  component,  although  the 
coefficient  would  be  difficult  to  determine.  Figures  2  and  3  have 
marked  square  law  behaviour  but  the  overall  behaviour  is  such  that 
one  could  not  trust  the  results. 


II.  FRICTION  AND  WINDAGE  BY  MECHANICAL  MEANS 


As  outlined  in  Chapter  III  Section  III,  simple  damped  harmonic 


motion  would  lead  to  an  approximate  period 

f.  Z«J~ 

‘/(wrJ-R') 

From  6  at  t  t  one  would  obtain 
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Since  the  amplitude  could  be  directly  measured  and  Gf  was  measured, 
all  values  necessary  could  be  determined. 


at 

o 

n 

ij 

the 

weight 

was 

17.5 

inches 

above 

datum 

at 

the 

weight 

was 

8.75 

inches 

above 

datum 

at 

t-00  , 

the 

weight 

wa  s 

3.25 

inches 

above 

datum 

The  radius  of  the  shaft  was  1.118  inches  (including  one  half  the 
thickness  of  the  wire).  The  linear  slope  of  the  figure  4  was  used 
to  estimate  the  frictional  force  and  in  turn  the  retarding  torque 
was  estimated  at  7.39  inch-pounds. 


Thus  the  net  torque  was  estimated  at  11.9  inch-pounds 
at  t-  o  : 

14.25  =  A,  e°  cos  (o') 

=  Ae 


The  period  was  observed  to  be  2.85  seconds  and  thus 
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-  zr/zz 

e  cos  ZTf 
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or  instead  of  inch- lb-sec  radian  and  inch- lb-sec  radian 
respectively,  J  =  0.223;  R=  0.149  foot- lbs- sec^  radian  and 
foot-lb-sec  radian  1  respectively. 
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TABLE  I 


MECHANICAL  LOSSES  DETERMINATION:  SUPPLY  FROM  D.C. 

D.C.  MACHINE  COUPLED  TO  A.C. 
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Power  Power  Constant  1.20  Amps. 
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APPENDIX  D 


ARTIFICIAL  DAMPING 


The  writer  altered  the  program  to  study  the  case  where  a 
torque  is  applied  tc  the  rotor  shaft  according  to  the  formula: 

Ta  =  -  k  (  e  -  e0 ) 

where  ©  is  the  actual  velocity 
e0  is  the  ideal  velocity. 


The  new  differential  equation  is 

4  Re  »  Le  -  K  (  @  -  ©a ) 

where  Te  is  the  developed  electromagnetic  torque  in  the  doubly 
excited  induction  machine. 


Then 
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With 

K  .  o-  9 65/4  R  _  0-149 

K4  R 

K  4  R  _  \<B‘  \&S  7  -  0-223 


O’ 

an  improvement  was  noted  but  stability  was  not  achieved.  For  the  above 
data  K  -  4.1252  (ft-lbs-sec  radian  )  which  is,  in  reality,  a  very 
small  value  in  comparison  to  what  can  be  achieved. 


It  should  be  noted  that  the  torque  given  by 

Ta  _  _*(©-©*) 

is  very  close  to  the  torque  developed  by  a  shunt  wound  direct  current 


.  v 


motor. 


The  accompanying  two  coding  sheets  demonstrate  the  changes 
required  to  study  the  phenomena  associated  with  the  externally  applied 
torque.  (See  statement  (15)  Appendix  B) . 
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APPENDIX  E 


CONSTANT  CURRENT  REGULATOR 


Shunt  field  heating  in  the  separately  excited  shunt  connected 
direct  current  motor  resulted  in  some  drift  in  shunt  field  current. 
Inductive  effects  and  possibly  other  unknown  or  undefined  functions 
such  as  armature  reaction  would  have  made  many  of  the  experiments 
invalid.  Accordingly  the  writer  sought  methods  to  regulate  the  shunt 
field  current. 


Assistant  Professor  K.G.  Clarke  kindly  loaned  the  writer  two 
power  transistors  and  one  zener  diode  and,  far  beyond  the  call  of 
duty,  designed  a  circuit  using  these  components  which  provided 
remarkably  constant  current.  The  circuit  as  designed  by  Mr.  Clarke 
is  included  herein  and  a  detailed  method  of  operating  the  device  is 
also  given. 


shunt 


field 


Fi 


9- 


I 


CyRfctNT  RmulATOR. 
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124. 


The  function  of  the  zener  diode  was  to  provide  a  comparatively 
small  reference  current  for  R2  and  the  function  of  T2  was  to  amplify 
this  reference  current  to  such  a  level  that  the  combined  current  would 
be  at  least  ten  times  the  base  current  of  Tt, 

With  a  drop  in  field  current  the  bias  across  the  transistor 
from  emitter  to  base  would  tend  to  decrease  and  the  mutual 
conductance  would  tend  to  rise  allowing  more  collector  current. 

Since  an  increase  of  R,  beyond  a  certain  maximum  would 
produce  an  emitter  potential  less  than  base  potential  and  result  in 
negative  base  current,  care  had  to  be  exercised  in  use  of  the  equipment. 
Consequently  the  following  procedure  was  adopted: 

(a)  The  value  of  R,  was  set  to  minimum. 

(b)  The  variable  voltage  supply  V,  was  applied  at  a  value 
which  could  not  have  delivered  more  than  desired  current 
to  the  field  windings. 

(c)  If  the  product  of  emitter-collector  voltage  and  collector 
current  exceeded  10  watts  V,  was  reduced,  if  not  10  watts 
the  value  of  R,  was  increased  until  the  product  was  10 
watts  or  the  emitter-collector  voltage  reached  20  volts 
whichever  occurred  first. 

It  was  observed  that  the  emitter  to  collector  voltage  could 
be  controlled  almost  exclusively  from  V,  while  R,  controlled  the  current 
almost  exclusively.  In  essence,  any  field  heating  effects  were 
cancelled  by  a  reduction  of  emitter-to-collector  voltage  and  this 
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voltage  had  to  be  kept  reasonably  high  in  order  to  avoid  much 
adjustment  of  the  apparatus.  So  long  as  the  emitter  was  higher 
in  potential  than  either  the  base  or  collector,  the  operator  could 
be  certain  that  the  field  current  was  constant. 
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APPENDIX  F 


SUMMARY  OF  ELECTRICAL  AND  MECHANICAL  DATA 

I.  ELECTRICAL  DATA 

With  the  rotor  removed  from  the  machine,  the  resistances 
and  reactances  of  each  phase  were  determined  as  follows: 

(a)  Three  single  phase  wattmeters  were  selected  having 
nearly  identical  impedances  in  order  to  avoid  unbalancing  the  three 
phase  system. 

(b)  The  star  point  of  the  rotor  was  made  accessible. 

(c)  All  three  wattmeters  had  their  potentials  referred 
to  the  star  point. 

(d)  From  data  collected  the  following  constants  were 
calculated : 
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The  voltage  was  kept  sufficiently  low  that  hysteresis  and 
eddy  current  losses  would  be  negligible,  the  phase  potentials  being 
a  little  under  four  per  cent  of  rated  rotor  phase  potentials. 

A  similar  procedure  was  used  to  determine  the  stator 
quantities.  In  this  instance  the  stator  phase  potentials  were  kept 
below  ten  per  cent  of  rated.  The  constants  calculated  for  the 


stator  are  as  follows; 
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The  writer  noted  that 
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A  variation  of  less  than  5  per  cent  is  observed.  The  coefficients 
of  self  inductance  are  proportional  to  the  square  of  the  turns  and 


the  voltages  are  proportional  to  the  turns.  Hence  the  voltages  are 

proportional  to  the  square  root  of  the  reactances.  If  the  stater 

voltage  is  220  volts  then  the  rotor  voltages  should  be  220  =  86.2 

'Jeso 

volts . 


The  writer  found  that  stator  potentials  of  224,  227  and  223  volts 
required  86.2,  83.1  and  86.8  volts  in  the  rotor  before  parallelling. 
This  would  indicate  an  error  of  perhaps  two  per  cent  which  is  entirely 
consistent  with  instrumental  error  and  personal  errors. 


II.  MECHANICAL  DATA 

With  the  procedure  outlined  in  Chapter  IV  Section  V  the 
following  constants  were  calculated: 

(1)  J  =  0.223  (F.P.S.  system) 

(2)  R  =  0.149  (F.P.S.  system) 
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Ill .  MUTUAL  INDUCTANCE 

The  mutual  inductance  was  determined  by  applying  the 

f ormu la 

The  above  formula  tends  to  neglect  leakage  which  would 
be  of  the  order  of  three  to  six  per  cent  and  m  could  be  in  error 
by  one  and  a  half  to  three  per  cent.  One  could  determine  the  mutual 
inductances  from  open  circuit  tests  of  the  rotor  but  it  has  already 
been  demonstrated  that  the  results  of  the  above  formula  ere  acceptable. 

The  value  of  m  calculated  was  0.04955  henries. 


IV.  SUMMARY 
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